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ABSTRACT 
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. Aims. We study the relations existing between fluxes emitted at CO (1-0) line, 60 and 100 /im wavelengths, B and soft X-ray 
wavebands for galaxies of all morphological types. The large set of data that we created allows to revisit some of known 
\£> relations existing between the different tracers of the Interstellar Medium (ISM): the link between the FIR flux and the CO 
f*"^ ■ line emission, the relation between X-ray emission in non active galaxies and the blue or FIR luminosity. 

Methods. Using catalogues of galaxies and works presented in the literature, we collected fluxes in FIR, 21 cm, CO(O-l) line 
and soft X-ray for two samples, consisting of normal and interacting galaxies respectively. Joining together these samples, we 
have data for a total of 2953 galaxies, not all observed in the four above wavebands. 
O Results. All the relations found are discussed in the frame of the star formation activity that is the link for most of them. 
£j We note that when an active star formation is present, it may link the galaxy fluxes at almost all wavelengths, from X to 
^2 ' microwaves. On the contrary, in early-type galaxies where the current star formation rate has faded out the X-FIR fluxes link 
disappears. This result obtained for early-type galaxies is discussed and explained in detail in the frame of a suitable theoretical 
model, obtained coupling chemo-dynamical N-body simulations with a dusty spectrophotometric code of population synthesis. 
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1. Introduction 

The observations of galaxies at various wavelengths, go- 
ing from radio to X-ray, allow to study the relation- 
ships existing between the various phases of the in- 
terstellar gas, and between gas, dust and stars. Some 
of these relations are already known since many years, 
such as th at between CO and far infrared (FIR) lu- 
minosities llSanders fc MirabeJll985l ISanders et aZ.lll986l 
ISolomon fc Sagell988l:lDevereux fc Young Ill99l|) . Others- 
connecte d with X-ray emission , have been studied more 
recently ^Griffiths fe Padovanil Il99(j | David et al 1 119921 
iRanalli et al\\200$ . 

At present, different tracers of the gas are known, such 
as millimetric lines for the cold molecular gas, the 21 
cm line for atomic hydrogen at ~100 K, IR bands for 
molecules at thousands of degrees, UV lines and X-ray 
emission for hotter gas. The dust distribution is traced also 
by FIR emission at 60 and 100 /im if the grains are warm 



Bregman_e^ 



Popescu et all 



1992) or at 170 /im, if they are colder 



2002). The diffusion of large archives of 
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observations at the above wavelengths (except for molec- 
ular lines) allowed in the last years the compilation of cat- 
alogues containing a huge number of galaxies. Using these 
catalogues and the works presented in the literature, we 
collected fluxes in FIR, 21 cm, CO (0-1) line and soft X- 
ray for two wide samples of normal llBettoni et a^jE0 03a') 
and interacting l|Casasola et a^Jl2004a|) galaxies. Joining 
together these samples, we have data for a total of 2953 
galaxies, not all observed in the four above wavebands. 
The fluxes measured with the different tracers allow now 
a study of the link existing between dust, gas and stars 
based on hundreds of galaxies. 

It is known that the fluxes emitted by a galaxy at very 
different wavelengths may be linked together by mean s 
of the star formatio n mechanism (see lDavid et all (^992); 
IRanalli et all l|2003h ). For instance, the formation of mas- 
sive stars generates the heating of the dust clouds in which 
they are embedded, by absorption of their UV radiation, 
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Fig. 1. The relation between the flux from the CO(l-O) line and that at 100 /im, plotted as a full line, with the 
representative points of the single galaxies plotted with different symbols. The crosses at the upper le ft of the figures 
repres ent the mean errors for the fluxes in the whole sample. For comparison the relation found by iBregman et all 
(1992) for a sample of early-type galaxies is shown as a dotted line. The early-type galaxies and early spirals, from E 
to Sb, are plotted in the left panel while the late spirals are plotted in the right panel. 



and produces a re-emission of this energy in the far in- 
frared. This process links the curren t star formation rate 
to th e IR emission at 60 and 100 fj,m IjThronson fc Telescol 
Il986|) . The ionizing radiation of stars may produce also the 
evaporation of the molecular clouds. Inside these clouds, 
where the particle density is great enough to produce 
a significant number of collisions between H2 and CO 
molecules, these latter are excited and produce photons, 
but in optically thick regions. The warming by the UV 
stellar light makes these regions less dense, making vis- 
ible the CO lines at their edge. Because of this mecha- 
nism, these lines are considered tracers of the cold molec- 
ular hydrogen that does not emit observable lines. The 
newly formed stars are also responsible of the X-ray emis- 
sion, produced by very massive stars, by core-collapse SN, 
and by high mass X-ray binaries. According to the above 
described mechanisms, we expect that galaxies with ac- 
tive star formation will have a far infrared emission, but 
also CO and X-ray emissions induced by the more massive 
stars, linked together by means of different relations. 

When the star formation decreases or vanishes, the far 
infrared emission decreases as well, but it may be fed by 
the stellar light absorbed and re-emitted in the infrared 
by dust (cirrus), while low-mass X-ray binaries and Type 
I SN contribute to the high energy galaxy spectrum. In 
addition, AGB stars, surrounded by dust, and the cool- 
ing flows of the interstellar medium ejected by supernovae 
may produce additional IR and X emission, between each 
other. 

To study the activity of the galaxies at different wave- 
bands, we collected data on galaxies starting from the orig- 



inal data of fluxes at 60, 100 /im, CO(l-O ) 2.6 mm and 
soft X - ray used to compile ou r catalogues l|Bettoni et all 
2003lj ICasasola et oiJl2004bft . The merging of the two 
above catalogues produces 1764 known values of far in- 
frared fluxes (1837 have 100 ^m flux), 391 soft X-ray 
fluxes and 434 values of the CO( T0) line luminosity . We 
extracted from LEDA catalogue ((Paturel et 997ft the 
values of the distance moduli, blue absolute magnitudes 
and morphological classification for all of them. 



Galaxies with evident sign of interactions or di s turbe d 
morphologies a c cording to the catalogues of lAxrJ lll966ft ; 
lArp &: Madord l)l987ft ; T^rontsov-Velvaminovl l)l959|) are 
1038. We shall refer to them as "perturbed galaxies". The 
remaining 1915 galaxies that appear neither morphologi- 
cally nor dynamically perturbed are called "normal galax- 
ies" . In our sample, we have 253 galaxies that have spec- 
tral classification of the nucleus and 231 of these appear to 
host an AGN (Seyfert 1, 2 or transition t ype, Seyfert 3 or 
Line rs) according to the classifica tions of iHo et all ( 1997ft 
and IVeron-Cettv & Veronl l|2003ft . Most part of the re- 
maining 2722 galaxies lacks of information about nuclear 
spectrum or have spectra of HII regions (22 starburst spec- 
tra). They are not included in any AGN catalogue and for 
this reason in the following discussion we refer to them as 
"non active galaxies" and to the others as "active galax- 
ies" . With all these data, we crossed the various tracers to 
understand and revisit the main relations existing between 
X, FIR, CO and B luminosities. 
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Fig. 2. Left: The X-ray luminosity plotted vs. blue luminosity, in solar units, for late-type, non-active galaxies. The 
relation corresponding to emission caused by discrete sources, indicated in equation (@J is plotted as a dotted line. 
Right: X-ray luminosity vs. FIR luminosity, in solar units for the late-type, non-active galaxies. The linear relation 
indicated in equation J7J) is plotted as a dashed line. 



2. Cold gas and warm dust 

The relations existing between different cold components 
of the ISM such as the molecul ar gas and the dust have 
been studied since many years llSanders fc Mirabe]|ll985t 
ISolomon fc Saeei ll98St iBreeman et al\\l99$i . Thev find 
that the global galaxy luminosity derived from CO (1-0) 
line is directly related with the flux at 100 /im. With our 
large sample we can now test these relations using galax- 
ies of different morphological types and activity or inter- 
action. 

In Figure H w e plotted the logarithm of the flux mea- 
sured from CO(l-O) line vs. the logarithm of the IRAS 
flux at 100 fim. In our plots, we have 193 galaxies with 
classification from E to Sb and 178 from Sbc to Sm. The 
relation found by iBregman et all (|l992^l for a sample of 
early-type galaxies, log Sco=log Sioo - 1-76, is also plot- 
ted as comparison, as a dotted line. 

The relations are evident, with this wider sample of 
galaxies. In these diagrams, active and non-active galax- 
ies appear mixed together without clear differences and 
have been plotted together. The same behaviour appears 
for interacting and non interacting galaxies, that are not 
distinguished in our plots. 

For all the galaxy types, we find: 



LogSco — 1-06 LogS 



too 



2.02 



(1) 



with a correlation coefficient of 0.74 and a r.m.s. of 0.37. 
In the above formula, Sioo is m m Jy and Sco is in Jy 
km/s. 

Similar relations exist between the CO fluxes and the 
FIR magnitudes, defined as: 

m FIR = -2.5 Log{2.58 S 60 + Sioo) + 22.25 (2) 



where Sgo and Sioo, the fluxes at 60 and 100 /zm respec- 
tively, are in mJy. We find for all the galaxy types: 



LogSco = 0.41 tufir + 6.86 



(3) 



with a correlation coefficient of 0.69 and a r.m.s of 0.40. 
The results are based on 179 early types and 170 late- type 
galaxies. For their similarity with Figures^these relations 
are not plotted in this paper. 

We note that irregular galaxies are not fitted by these 
relations but have a wide spread. In our sample there are 
just 10 galaxies and their representative points have been 
not plotted in Figure 



3. X-ray component. 

We are interested to understand what relations exist be- 
tween Lx, the X-ray luminosity, and the other global 
galaxy properties. From the literature, it is known the 
existence of a proportionality between Lx produced by 
discrete sources and Lb, the blue luminosity of the whole 
galaxy . This relation h as been studied by ICiotti et 



||1991I) and compared bv lBeuing et all l)1999t) with soft X- 
ray fluxes measured by ROSAT satellite. It appears that 
late-type galaxies have a global X-ray luminosity directly 
proportional to Lb, while early- type systems are domi- 
nated by emission produced by hot diffuse gas and their 
Lx is proportional to the square powe r of th e blue lu- 
minosity, as discussed by iBeuing et all l|l999h . For this 
reason, the early and late-type galaxies are discussed sep- 
arately. 




Fig. 3. Left: The X-ray luminosity plotted vs. blue luminosity, for early-type, non-active galaxies. The relation cor- 
responding to emission caused by diffuse gas, valid for early-type galaxies and indicated in equation is plotted 
as a full line. Right: X-ray luminosity vs. FIR luminosity for the early-type galaxies. The linear relation for late-type, 
non-active galaxies indicated in equation (JJJ is plotted as a dashed line. 



3.1. Late- typ e ga laxies 

With our data, the X-ray luminosity of galaxies with mor- 
phological type later than Sb can be fitted by a linear re- 
lation as a function of Lb (dotted line in Fig.|2| left panel). 
The direct proportionality is expressed by the equation: 



LogLx = LogLs — 3.85 



(4) 



with a r.m.s. from observed data of er=0.61 based on 63 
galaxies. In this formula and in the following, all the lu- 
minosities are expressed in solar units. 

If, instead of the blue luminosity, we use the galaxy 
area Df. pc , calculated from the apparent diameter mea- 
sured at the 25 mag arcsec 2 isophote and converted in 



kpc 2 , we discover that the relation is still present, 
with a larger spread. It becomes: 



LogLx = LogD 



2 

kpc 



3.83 



but 



(5) 



(cr=0.80) for a sample of 64 gala xies. 

A relation similar to that of ICiotti et 



[ 199 3 has 

been found by some authors ^Griffiths fc Padovanii ri990: 
iDavid et flj]ll992t iRanalli et all 1200.31). but usinTfiO M m 
fluxes or FIR luminosities. The values of Lfir. are calcu- 
lated using the formula: 



to Lx oc Lpj R found bv lDavid et all Jl992) using fluxes 
between 0.5 and 4.5 keV. Forcing the relation to a linear 
proportionality between Lx and Lfir we find: 



LogLx = Log L fir - 3.18 



(7) 



with a a of 0.47, based on 147 galaxies. This relation is 
plotted as a dashed line in the right panels of Figures from 

EltoH 

We note that the B and FIR luminosities are also con- 
nected in late-type galaxies by means of a linear relation 
fitted by: 

LogLpiR. — LogLs — 0.38 (8) 

with a r.m.s. =0.5. This equation, inserted into the relation 
(0J gives: 

LogLx = LogLpiR - 3.47 (9) 

similar to the result of equation Q and to that found by 
IRanalli et all l|2003h . This is an independent way to con- 
firm our results and to verify the existence of a global link 
between Lpm, B light and X-ray emission. The connec- 
tion between B luminosity or galaxy area and X or FIR 
luminosities will be discussed in Section [SJ 



LogL fir = 2.59 + Log(2.58 S 6 



2Log d (6) 3.2. Early-type galaxies 



where Lrir is in solar luminosities, fluxes are in mJy and 
the galaxy distance d is in Mpc. 

From our data it is possible to find a relation between 
Lx and Lrir that fits the values of late-type galaxies. We 
found Ly (x Lp?p similar to the Lx oc L FIR found by 
IRanalli et all l)2003|) for fluxes between 0.5 and 2 keV and 



When the early-type galaxies are considered in the above 
described relations involving X-ray emission, the correla- 
tions become less evident. Considering soft X-ray and B 
luminosities, we find a relation: 



LogLx = 2 LogLR — 13.57 



(10) 




Fig. 4. Left: The X-ray luminosity plotted vs. blue luminosity, for active galaxies, separated according the type of 
activity. The relation corresponding to emission caused by diffuse gas, valid for early-type galaxies and indicated in 
equation i|l(Jfl is plotted as a full line. Right: X-ray luminosity vs. FIR luminosity for the active galaxies. The linear 
relation for late-type galaxies indicated in equation (JJJ is plotted as a dashed line. 



(cr=0.73) based on 224 galaxies and plotted as full line in 
Fig. 01 left panel. The above formula agrees with the ex- 
pected relation for X-ray emission comin g from hot diffuse 
gas, as discussed bv lBeuing et al\ l|l999F) . 

The relation still hold if Dl pc (kpc 2 ) is used. It be- 
comes: 

LogL x = 2LogD 2 kpc + 1.51 (11) 

(cr=0.85) for 226 early-type galaxies from E to Sb. 

Many galaxies with high blue luminosity, indication of 
high masses and of a recent star formation, lie quite far 
from the mean line, with a behaviour different than that 
of late-type galaxies. 

If the X-ray fluxes are compared with FIR luminosity, 
the disagreement with the behaviour found in late-type 
galaxies is more evident. The plot Lx vs. Ijfir for early- 
type galaxies shows the representative points of the galax- 
ies above the relation Q for late-type galaxies (Fig|3l 
right panel). To understand this apparent disagreement, 
we should use a theoretical analysis of the far infrared 
emission, as explained in the next Section 4. 

3.3. Active galaxies 

Active galaxies (Seyfert 1, Seyfert 2 and Liners) have X- 
ray, B and FIR fluxes that are not linked together. This 
happens because, to the emission mechanisms stimulating 
the light emission at the different wavebands described for 
non active galaxies, adds an X-ray emission coming from 
nucleus. In fact, the points representative of these active 
galaxies are spread in the plot over the discrete sources 
line and around the diffuse gas line (see Fig0J left side). 
In the Lx-Lpm diagram (FigfJ] right side) the spread 



is similar to that of early-type galaxies plotted in FigOU 
but we separately plotted the active galaxies because of 
the particular nature of their X-ray emission, due to the 
nuclear contribution. 



4. Modelling Lx, Lb and Lfir of early-type 
galaxies 

To cast light on the nature of the relations observed be- 
tween Lx, Lg and Lpm for early-type galaxies, one has 
to consider the various components of a galaxy (stars, gas 
and dust) and to understand their mutual interactions as 
far as the spectral energy distribution (SED) is concerned. 
There are two basic schemes to model the formation and 
evolution of early type galaxies: (1) the semi-analytical 
models on which a great deal of our understanding of the 
chemo-spectro-photometric properties is derived, and (2) 
the N-Body Tree-SPH simulations which, in contrast, have 
been only occasionally used to study spectro-photometric 
properties of early type galaxies. In the following part of 
this section we will proceed as follows. First we will anal- 
yse the drawbacks of semi-analytical models, in particu- 
lar dealing with the calculation of the infrared emission of 
early-type galaxies. Second, we will discuss how dynamical 
simulations and a dusty spectrophotomctric code, when 
mixed together allow to move a step forward in the cal- 
culations of the SEDs properties. Third, we will show in 
detail how our model has been built and the coupling be- 
tween dynamics and dusty population synthesis has been 
done. 



6 



G. Galletta et al.: ISM tracers 



4.1. The semi-analytical models and their drawbacks 

The semi-analytical models approximate a galaxy to a 
point mass system in which gas is turned into stars by 
means of suitable recipes for star formation and heavy 
elements are produced by stellar nucleosynthesis and stel- 
lar winds/explosions. The standard evolutionary popula- 
tion synthesis technique (EPS) is usually applied to de- 
rive the SED of the galaxy, with models able to explain 
many global features of e arly type galaxies, as amply de- 
scribed by many au t hors JArimoto fc Yoshii 111 987. 1989; 
Bressan et all Il994t iGibsonl Il997t iTantalo et all Il996 , 
1998|) . There are three important and problematic issues 
of these models to be discussed for our purposes. 

First, to determine the age at which the gal actic 
wind sets l|Larson Ill974t lLarson fc Dinerstenlll975|) . we 
need some hypothesis about Dark and Baryonic Matter 
with their relative distributions, and about the heat- 
ing and cooling efficiency of the various mechanisms, 
to properly evaluate the total gravitational potential 
well and to describe the thermal history of the gas. 
In this scheme it comes out that the galactic wind oc- 
curs typically for ages taw < 1 Gyr, later in a mas- 
sive early- t ype galaxy and much earlier in galaxies of 
lower ma,ss JArimoto fc Yoshii 111 987lll 989l iBressa.n et al 
1994 [lGibsoJl997MTantak> et aZJIl996lll998tlChiosi et al 
1998). The maximum duration of star forming activity 
follows therefore in these models the trend At$F Mq- 
This trend of the SFH is, however, contrary to what 
required by the observes trend of the a-enhancement 
for early type galaxies, which implies that the maxi- 
mum duration of the star forming activity should de- 
crea se when the galaxy mass increases f Ai.sF oc M^ 1 ) 
(seelBressan et allll99(l iKuntschner Il2000t iTrager et all 
l2000alhl: ITantalo fc Chiosi ll2004HThomas et allhOOFl for 
more details on the enhancement in a-elements and the 
SFH of early- type galaxies). 

Second, after the galactic wind phase, star formation 
does no longer occur and the evolution is merely passive. 
However, AGB and RGB stars continue to loose gas in 
amou nts that are c omparable to those before the galactic 
wind l)Chiosi l2000|) . What is the fate of this gas? One may 
imagine that the large amount of gas lost by stars will ex- 
pand into the Dark Matter halo and heat up to an energy 
overwhelming the gravitational potential, it will escape 
the galaxy. Most likely a sort of dynamical equilibrium is 
reached in which gas is continuously ejected by stars and 
lost by the galaxy. It may happen therefore that some 
amount of gas is always present in the galaxy. The ques- 
tion is not trivial because if an early type galaxy is free of 
gas and contains only stars, the SED is expected to drop 
off long-ward of about 2/im and no IR emission should 
be detected. Howeve r , as a lr eady pointed out lo ng ago by 
iGuhathakurta et ail <)l986(klKnapp et al\ l|l989f) (see also 
Fig. El, many early-type galaxies of the local universe emit 
in the IR. The origin of this flux in the MIR/FIR is likely 
due to dust present in a diffuse ISM which, heated up by 
the galactic radiation field, emits at those wavelengths. 



Therefore to match the IR emission one has to allow for 
some amount of diffuse ISM. An interesting question to 
rise is therefore: how much gas can be present today in 
an elliptical galaxy and how is it distributed across the 
galaxy? Even if we can correctly estimate the amount of 
gas ejected by stars, the fate of this gas goes beyond the 
possibilities of classical semi- analytical models. 

As a third point, note that when we fold many SSPs to 
calculate a galaxy SED using the classical EPS technique 
we simply convolve their fluxes with the SFH of the galaxy. 
Many classical spectrophotometric semi-analytical models 
of galaxies are built in this way: there is no dust at the level 
of SSPs a nd again no dust at the leve l of the galaxy model 
(see e.g. lArimoto fc Yoshii I Il987t lArimoto fc Tarrab I 
199(i iBruzual fc CharloTl Il99.lt ITantalo et all Il99fit 
Kodama fc Arimoto Ill997t ITantalo et oi.lll998t iBuzzoni I 
20021 12005|) . To calculate the emission by dust, a higher 
level of sophistication of the model is required. Indeed one 
has to develop a model in which the sources of radiation 
and the emitting/absorbing medium are distributed, to 
face and solve the problem of the radiative transfer simu- 
lating in a realistic way the interactions among the various 
physical components of a galaxy. Among recent models of 
this k ind a re those bvlSilva et all l|l998l) . lDevriendt et all 
!|1 999l) and lTakagi e t all l|200,^ . 

4.2. Improving upon semi-analytical models 

Two drawbacks of the semi-analytical models concern 
therefore: (1) the description of galactic wind, which is 
supposed to occur within a finite time interval and (2) 
the star formation history that is reversed allowing longer 
SFH for more massive galaxies. These two problems, com- 
bined with a lack of geometrical information about the 
distribution of gas and dust, make semi-analytical mod- 
els not suitable to calculate properly the IR emission of 
early type galaxies. To improve upon them we need to 
use the results obtained from dynamical simulations. They 
have shown to be able to properly model the ejection of 
gas by the galaxy as a sort of continuous process, tak- 
ing place whenever a gas particle heated up by various 
mechanism has acqui red a velocity greater than the es- 
cape velocity (see e.g.lCarraro et g,/Jll998t Fkawata iboOlt 
Snringel H200lilGhiosi fc Carraro 112 0021. They are able to 
reproduce the SF history of early-type galaxie s both in 
the context of the monolithic collapse scenario l|Kawata I 
l200lHChiosi fc Garraro l|2002l) and recently in the context 
of hierarchical scenario l)De Lucia et all\200(^ . 

Finally, the galaxy is no more a mass point, but a fully 
three-dimensional structure of the galaxy is available with 
spatial distribution of stars and gas. 

iMerlin fc Chiosi 1 l)2006|) . with the aid of N-Body Tree- 
SPH simulations based on quasi-cosmological initial con- 
ditions in the standard-Cold Dark Matter scenario (S- 
CDM), modelled the formation and evolution of two early- 
type galaxies of different total mass (Dark + Baryonic 
Matter in the cosmological proportions 9:1). The total 
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Tab le 1. Initial paramete rs for the dynamical simulations 
of lMerlin fc Chiosi 1 1)2006|) in the Standard CMD scenario. 
Masses are in units of 10 12 M Q , radii are in kpc and ages 
are in Gyr. 



Mo del 


A 


T"> 

D 


Cosmological background 


o-ODM 


b-OJJM 


Initial redshift 


50 


53 




1 


i 


Nq — 5U kmMpc s 


50 


50 


Gas particles 


13719 


13904 


IjlJM particles 


1 QC Q ^ 

lODOO 


L6l 10 


rp„i 1 TV T 

lotal Mass 


1.62 


0.03 


Initial baryonic mass fraction 


0.10 


0.10 


Present gas mass 


0.062 


0.0004 


Present star mass 


0.091 


0.0029 


^Istar 1 ^ibaryons 


0.556 


0.82 


Initial radius 


33 


9 


Half-Mass radius of stars 


7 


1 


Half-Mass radius of DM 


52 


15 


Effective radius of stars 


5.2 


0.8 


Present virial radius 


300 


41 


Axial ratio b/a (stars) 


1.08 


1.04 


Axial ratio c/a (stars) 


1.07 


1.00 


Axial ratio b/a (Dark Matter) 


1.14 


1.14 


Axial ratio c/a (Dark Matter) 


1.17 


0.96 


Age of the last computed model 


13 


5 



masses under considerations are 1.62 x 1O 12 M0 (Model 
A) and 0.03 x 10 12 M o (Model B). The galaxies have been 
followed from their separation from the global expansion 
of the universe to their collapse to virialized structures, 
the formation of stars and subsequent nearly passive evo- 
lution. They are followed for a long period of time, i.e. 13 
Gyr (Model A) and 5 Gyr (Model B). In any case, well 
beyond the stages of active star formation which occurs 
within the first 3 to 4 Gyr (see below). The models take 
into account radiative cooling by several processes, heat- 
ing by energy feed back from supernova explosions (both 
Type I and II) and chemical enrichment. All the models 
conform to the so-called revised monolithic scheme, be- 
cause mergers of substructures have occurred very early 
in the galaxy life. Some parameters and results of the two 
models are summarized in Tabled Note that the shape 
of the resulting galaxies is nearly spherical both in Dark 
Matter and stars. 

The third drawback of classical semi-analytical model 
was the lack of the description of the dusty component, 
that for our purposes needs to be included. The semi- 
analytical chemo-spe ctro-photometric model developed by 
Pio van et Mil (|2006b) allows us to overcome this issue. It 
takes into account not only the geometrical structure of 
galaxies of different morphological type, but also the ef- 
fect of dust in converting t he UV and Optical light in 
far IR radiation. In brief the lPiovan et all l)2006bh model 
follows the infall scheme, allows for the onset of galactic 
winds, and contains three main components: (i) the dif- 
fuse interstellar medium composed of gas and dust whose 
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Fig. 5. Top panel: The masses of stars (continuous line) 
and gas (dashed line) for the dynamical model of 1.6 • 
1O U M baryonic mass as a function of the number of 
virial radii. The virial radius for this model is about 300 
Kpc. Bottom panel: The same of the upper panel, but for 
the dynamical model of baryonic mass 3.5 ■ 1O 9 M . For 
this model the virial radius is about 42 Kpc. 



emission and extinction prop erties have been studied in 
detail bv lPiovan et al\ l)2006a|) . (ii) the large complexes of 
molecular clouds in which new stars are formed and (iii) 
the stars of any age and chemical composition. The to- 
tal gas and star mass provided by the chemical model are 
distributed over the whole volume by means of suitable 
density profiles, one for each component and depending 
on the galaxy type (spheroidal, disk and disk plus bulge). 
The galaxy is then splitted in suitable volume elements to 
each of which the appropriate amounts of stars, molecular 
clouds and interstellar medium are assigned. Each elemen- 
tal volume absorbs radiation from all other volumes and 
from the interstellar medium in between. The elemental 
volume also re-emits the absorbed light and produces ra- 
diation by the stars that it contains. On the other hand, 
the star formation , the initial m a ss func tion, the chemical 
enrichment of the Piovan et model are much 

similar to t h ose bv iBres san et all l)l994|) ; iTantalo et all 
<|l99dll998jklPortinari et all (|l998|) . 

4.3. Coupling dynamical simulations and dusty 
population synthesis models 

The description of an early-type galaxy as far as predicting 
its spectro-photometric infrared properties can be there- 
fore realized with a suitable combination of dynamical and 
spectro-photometric approaches. Coupling the dynami- 
cal models with spectro-photometric synthesis requires a 
number of steps that deserve some remarks. 
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Fig. 6. The averaged density profiles of gas and stars for 
the model of 1.6 • 1O 11 M0 baryonic mass at the age of 13 
Gyr. The continuous lines are the best fit obtained. 



4.3.1. Radial density profiles. 

Fig. [3] shows the cumulative distribution of gas and stars 
as a function of the radial galactocentric distance normal- 
ized to the virial radius for model A (top panel) and model 
B (bottom panel). The gas is generally distributed in the 
external regions of the galaxy and steeply decreases in- 
ward. In contrast the stars are more concentrated toward 
the centre. The gradients in the spherically averaged star- 
and gas- content provided by the dynamical models are the 
primary info rmation t o load into the spectro-photometric 
code of IPiovan et all l|2006tJ) . They allow us to infer the 
amount of gas contained within a given radius or within 
a given aperture. We fix the total dimension of portion 
of the average model producing the IR flux at a diameter 
Dg a i = 25 kpc, consistent with the mean galaxy size of 
the observed sample. 

As the spectro-photometric code of IPiovan et all 
suited to describe early-type galaxies is written 
in spherical symmetry, we have to derive suitable spher- 
ical distributions for the density of stars and gas to be 
used into the model. The task is facilitated by the nearly 
spherical shape of the dynamical models. To this aim, we 
consider the sphere of radius R ga i centred at the centre of 
mass of the stellar component. The sphere is then divided 
in a number of thin spherical shells whose derived average 
density of stars and gas is shown in Fig. |SJ Even if the 
centre of mass of the star and gas distributions may not 
be exactly coincident, this not relevant here, so that the 
same coordinate centre can be used for both components. 

In order to secure a smooth behaviour at the galaxy ra- 
dius Rg a i the star and gas density profiles are represented 
by the law: 



PQi 



1+ -T 



21 ~1* 



(12) 



where "i" stands for "stars" or "gas", r\ are the cor- 
responding core radii. The above representation is more 
suited to our aims than the classical King law. The fits 
are shown in Fig. (solid lines) . They are normalized in 
such a way that the integral over the galaxy volume cor- 
responds to the amount of gas contained inside R ga i- 

4.3.2. Star formation rate. 

In the dynamical models, the period of intense star forma- 
tion, during which most of the star mass is built up, is con- 
fined within the first 3 to 4 Gyr. In Model A this is followed 
by a long tail of minimal stellar activity which continues 
forever. If this activity would be real, we would expect a 
background of young stars giving rise to a significant emis- 
sion in the UV-optical region up to the present, which is 
not compatible with the observed spectr a of typical early- 
type g alaxies. As already pointed out bv lMerlin fc Ghiosi I 
(2006J) this minimal stellar activity is an artefact of the 
poor mass-resolution for the baryonic component, in other 
words the low number of particles considered in the nu- 
merical simulations. To cope with this, we simply set to 
zero the star formation rate when only one or two star 
particles are involved. This is equivalent to cut the star 
formation rate for ages older than about 5 Gyr. The prob- 
lem does not occur with model B simply because the last 
computed models is at 6 Gyr. 



4.3.3. Checking dynamical models against 
chemo-spectro-photometric models. 

To this aim we plug the star formation history 
(SFH) of dynamical models into the chemical code of 
IPortinari et all £l998). The closed-box approximation is 
adopted. The total baryonic mass of the chemical models 
is the same as in the dynamical ones. Equally for the initial 
mass function of the stars composing each star particle: 
iKrouoa 1 l|l998l) in our case. In Fig. we show the results 
obtained by inserting the SFH of Model A into a classi- 
cal chemical model with total baryonic mass Mb equal 
to 1.6 • 10 11 M Q . The top panels display the adopted SFH 
(left) and the gas metallicity of the chemical model, re- 
spectively. The bottom left panel shows the temporal vari- 
ation of the star mass M sta r and gas mass M gas , whereas 
the bottom right panel shows the ratios M star /Ms and 
M gas /Mb for both the dynamical (thin lines) and chemi- 
cal model (thick lines). The agreement is very good thus 
confirming the internal consistency between the descrip- 
tions of the same object. We also show the amount of 
gas at 13 Gyr contained in the whole galaxy for both the 
dynamical (heavy dots) and the classical chemical mod- 
els (open circles) and the amount of gas contained inside 
Rgai (open squares). Indeed there is little gas left over 
inside the 25 kpc radius region. Similarly in Fig. |H1 we 
show the results obtained inserting the SFH of Model B 
into a classical chemical model with total baryonic mass 
of 3.5 • 10 9 M Q . The only difference is that the maximum 
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third, that we can safely apply t he population synthesis 
technique of lPiovan et all ( 2006bh . 



Fig. 7. Basic quantities of the chemical models for a pro- 
totype early-type galaxy of 1.6 • 1O 11 M baryonic mass as 
function of the age. Top-left panel: The star formation rate 
obtained from the chemo-dynamical model of 1.6 • 1O 11 M 
baryonic mass (SFR). Top-right panel: The evolution of 
the metallicity Z . Bottom-left panel: The masses for both 
the classical and the chemo-dynamical model of baryonic 
matter (thick and thin dashed lines), stars (thick and thin 
continuous lines) and gas (thick and thin dotted lines). 
The mass of gas and stars inside 25 Kpc is also shown 
(open circles) for the age of 13 Gyr. Bottom-right panel: 
The fractional masses for both the models of stars (thick 
and thin continuous lines) and gas (thick and thin dotted 
lines). 
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Fig. 8. The same as in Fig.[7|but for the model with total 
baryonic mass of 3.5 ■ 10 9 M Q . 

age of the dynamical model is 5 Gyr. This cross checking 
of the models is particularly significant because: first, it 
secures that the results of the analytical models fairly re- 
produce those of the dynamical simulations as far as some 
important features are concerned; second, it secures that 
we can safely use the result of chemical models to prolong 
the evolutionary history of Model B up to the present; 



Knowing the amount of gas, we need to specify the 
fraction of it in form of dust to finally be able to derive 
the whole SED from X to FIR and look for relationships 
between the luminosity in the X, B and FIR pass-bands we 
want to interpret. Our models, both semi-analytical and 
chemo-dynamical, are not suitable to describe the evolu- 
tion of the compositions and abundances of both gas and 
dust phases. The relative proportions of the various com- 
ponents of the dust would require the detailed study of 
the evolution of the dusty environment and the co mplete 
information on the dust yields, as in the models of lDwekl 
lll998tl2005h . This would lead to a better and more physi- 
cally sounded correlation between the composition of dust 
and the star formation and chemical enrichment history 
of the galaxy itself, however at the price of increasing the 
complexity and the uncertainty of the problem. 

The key parameter to calculate the amount of dust is 
the dust-to-gas ratio, defined as 5 — Md/Ma, where Md 
and Mh are the total dust and hydrogen mass, respec- 
tively. For the Milky Way and the galaxies of the Local 
Group, 8 is estimated to vary from about 1/100 to 1/500 
and typical values 6 = 0.01, 8 = 0.00288 and 8 = 0.00184 
are used for the Milky Way (MW) and the Large and 
Small Magellanic Clouds (LMC and SMC). These dust- 
to-gas mass ratios describe a decreasing sequence, going 
from the MW to the LMC and SMC. Since these galax- 
ies also describe a sequence of decreasing metallicity, a 
simple assumption is to hypothesize 8 oc Z in such a 
way to match the approximate results for MW, LMC and 
SMC: 8 = Sq (Z/Zq). This relation simply implies that 
the higher is the metal content of a galaxy, the higher is 
the abundance of grains per H atom. However, the metal- 
licity difference does not only imply a difference in the ab- 
solute abundance of heavy elements in the dust, but also a 
difference in the comp osition pattern as a function of the 
star formatio n history iDwek I jl998l,E005|) . Despite these 
uncertainties ( Devriendt et aiJll999h . the relation 8 ex Z 
is often ado pted to evaluate t he amount of dust in galaxy 
models (e.g. ISilva et qjJll998l) by simply scaling the dust 
content adopted for the ISM of the MW to the metallicity 
under consideration. 

The 1.6 • 10 n M© and 3.5 • 10 9 M Q galaxy models reach 
an average metallicity of solar and slightly more than twice 
solar, respecti vely. To describe them w e have adopted the 
description of IPiovan ~ ll2006albh where a model of 
dusty ISM taking into account different metallicities is 
built. The problem however remained unsettled for metal- 
licities higher than the solar one, where relative propor- 
tions holding good for the MW average diffuse ISM model 
have been adopted and the amount of dust scaled with 
8 oc Z. Therefore, for the 1.6 • 10 11 Mq galaxy with solar 
metallicity the MW diffuse ISM model has been adopted 
(8 = 8q), while for the 3.5 • 10 9 M Q model we followed the 
8 oc Z relation, using the MW average pattern of dust 
composition. 
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The connection between the results of this model and 
the observed diagrams are discussed in the following sec- 
tion. 

5. Discussion 

Our data confirm and extend the previous relations ex- 
isting between various tracers of the ISM in galaxies of 
different morphological types. 

In the literature the relation found by iBreeman et all 
(1992) between Sco and Sioo indicates a direct propor- 
tionalit y (slope=l) betwe e n the two fluxes and differs from 
that of ISolomon fc Sage] l)l988j) . that exhibits a steeper 
gradient. Our relation Q agrees quite well with the pro- 
portionality found bv lBregman et al . (1992). the slope we 
found being equal to 1.06. The similarity between the 
two curves in Figure Q is evident. This link derives, as 
described in the introduction, from the excitation of gas 
clouds by the currently forming stars and by the warming 
of the dust present in the galaxy. 

5.1. Late-type galaxies 

In late- type galaxies (t>Sb) our data show the existence of 
a linear relation between soft X-ray fluxes and other indi- 
cators of recent and current star formation, such as the B 
and FIR luminosity respectively (equations 0] and UJ. This 
is known since the firs t X-ray observations of large samples 
iFabbiano efoll {l992) and this connection between B and 
X-ray luminosity in late type galaxies has been interpreted 
as due to the contribution of discrete X-ray sources, whose 
numb er is proportional to the quantity of al ready formed 
stars JCiotti et aZJIl99ll iBemng et al\ \l§§§\. The recent 
work of iKim fc Fabbiano I 1)2004) . that is able to resolve 
the single X-ray binaries in 14 galaxies, indicates that the 
X-ray luminosity produced by discrete sources is related 
to B luminosity by a similar relation, with an intercept 
value of -3.63, similar to our -3.85 of equation^] 

In addition to the interstellar radiation, that is propor- 
tional to the number of already formed stars, the X-ray 
emission is produced also by HII regions, where ther e is an 
ongoing vigorous star formation l)David et a/Jll992h . This 
latter contribution appears more evident in FIR light and 
may explain the existence of a similar linear relation be- 
tween Lx and Lfir- 

5.2. Early-type galaxies 

In early type galaxies the behaviour of these relations 
is quite different. For most of these galaxies, the star 
formation is exhausted and it may be present in a few 
of them, eventually fed by gas accretion phenomena. 
Different mechanisms have been suggested to explain the 
X-ray emission in this kind of galaxies. In particular the 
main ones are the thermal emission due to hot ISM and 
the emission generated by a relatively old population of 
end objects of stellar evolution, composed by Type I su- 
pernovae remnants and low-mass X-ray binaries not yet 
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Fig. 9. Lfir vs. Lx diagram as in Figure 3, right panel, 
but with the values deduced for the two template models 
of early-type galaxy (filled circles). The diamonds repre- 
sents the E — SO galaxies, while the triangles the SOa — Sb 
galaxies. The masses of the two theoretical models are also 
indicated. 

evolved. In particular, for the fainter galaxies the X-ray 
emission is compatible with discrete sources and seems 
to be dominated by compact accreting systems, while for 
the brighter objects the emission from hot diffuse gas still 
present in the gal actic potential well is present as addi- 
tional component i Beuing et a/ J Il99flh . The number size 
of this population of relatively old objects is well repre- 
sented by the total blue luminosity of the galaxy. For this 
reason the X-ray fluxes are still linked in early type galax- 
ies with the total blue luminosity, representing the more 
recent part of the history of star formation in the galaxy. 

In the FIR however, since the star formation in most of 
these systems is almost exhausted, mechanisms different 
from the emission from warm dust heated by the newly 
born stars predominate. The FIR emission comes from 
circumstellar dusty shells around AGB stars and from an 
interstellar medium due to the outflow of dusty gas from 
AGB and RGB stars, as it has been described in Sect.|U 

The key point to interpret the observed trends is that 
we deal with an emission coming from a more or less small 
amount of dust distributed over all the galaxy and heated 
by an average interstellar radiation field due to all the stars 
of any age. The situation is quite different from what hap- 
pens for instance in starburst galaxies where high optical 
depth dusty regions reprocess the light coming from newly 
born stars embedded in the parental environment. We can 
therefore conclude that in most of our early-type galax- 
ies the mechanism of IR emission is not strictly related 
to the star formation and the link between the younger 
generations of stars and dust emission is lost. For these 
reasons one may expect that the soft X-ray luminosity in 
early type galaxies is traced by the total blue luminosity 
but not by the FIR luminosity. With the end of the star 
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formation, the far infrared emission of these galaxies has 
faded out and an early type galaxy with the same Lx of 
a late type will have a lower LpjR. This could explain the 
location of the points in Fig. [3] (right panel) , on the left 
side of the linear relation. 

To check if this interpretation is correct we try to apply 
the detailed chemo-dynamical spectrophotometric model 
described in the previous section, in such a way to estimate 
the luminosities produced by the stars in connection with 
the various phenomena present inside the galaxy, taking 
into account the contribution by dust as well. Since the 
theoretical model can not derive the Lx luminosity, we 
proceed in the following way. 

The luminosities Lb and Lfir are directly derived 
from the model. Then, we assume that the X-ray produc- 
tion of these galaxies is proportional to Lb according to 
our relation (|1C)|> . In this way we may estimate the ex- 
pected X-ray flux and define a representative point in the 
Lfir. vs Lx plot. 

We start considering two template models, in which 
all the parameters are fixed using the clues coming f rom 
the dynamical simulations of iMerlin fc Chio"sTI l|200fJ) . as 
described in Sect.0] The King profiles represented in Fig. 
|H|are similar for all the components, with -7 stars — Igas — 
1.5 and rf ars ~ r 9 c as ~ 0.5 Kpc, while the dimension 
of the galaxy is an average one corresponding to most of 
the galaxies available in the catalogue. The SFH is ex- 
actly the one obtained by the dynamical simulations. The 
two values of Lfir and Lx obtained for the 3.5 • 10 9 M Q 
and 1.6 • lO n M0 baryonic mass models are plotted in Fig. 
|5| The more massive galaxy fits well into the region de- 
fined by the really observed galaxies, while we can no- 
tice as the model of smaller mass, even if falling up to 
the linear relationship as we could expect, belongs to a 
region not covered by the observed data. The calculated 
levels of emission Lx and L fir of this galaxy are very low 
and for this reason they belong to a region where we do 
not have enough observations. The weak Lfir emission of 
this galaxy can be explained by the dynamical evolution 
in which almost all the gas is consume d to form stars and 
the g alactic winds are very efficient Csee lChiosi fc Carraro I 
120021 for more details about galactic wind in low mass 
galaxies). Therefore, even if the trend of this galaxy is 
the expected one for early-type galaxies (the model stays 
above the linear relation), nothing safer can be said, be- 
cause we lack observed data in that region of the diagram. 

Much more interesting is the model of higher mass. 
The calculated luminosities of the model, with its exhaus- 
tion of the star formation, seem to agree well with the 
observations of early-type galaxies. However, the model 
needs to be checked against other possibilities, at the pur- 
pose to understand the way in which the various param- 
eters of the model influence the spreading of early-type 
galaxies into the observational data. First of all we have 
to check the effect of the geometrical parameters and of 
the masses of stars/gas. 
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Fig. 10. Models of the 1.6 • 10 11 M Q baryonic mass galaxy 
in the Lpm vs. Lx diagram at varying the galactic radius 
(filled circles), the mass of stars (filled diamonds) and, 
finally, the mass of gas (empty circles). We represented: 
E galaxies (empty diamonds), E/SO galaxies (pluses), SO 
galaxies (crosses), SOs galaxies (empty squares), Sa — Sb 
(empty triangles). 

5.2.1. The galactic radius 

In fig.EDwe show the model of 1.6TO 11 M baryonic mass 
at varying the galactic radius, keeping the galactic center 
in the center of mass of the stellar component. The radii 
taken into account range from 6 Kpc to 50 kpc. All the 
other parameters are fixed. Four models are represented 
(filled circles) and connected by a continuous line and the 
smaller and bigger models are marked using an arrow. For 
larger radii we observe an increase of both L fir and Lx , 
with a more emphasized increase in Lx ■ Since the density 
profile is unchanged, both the increases in luminosity are 
simply due to the bigger amount of material considered 
taking into account larger radii in the dynamical simula- 
tion. The stronger increase in Lx than Lfir can be sim- 
ply explained. Lx is linearly related to Lb, that is directly 
connected to the stellar luminosity. The stellar component 
is more massive and more concentrated toward the cen- 
tre than the gaseous one (Fig.0 upper panel). It follows 
that at increasing radius we introduce into the models 
more stars and more gas, but the added amount of stars 
is bigger than the gaseous one, shifting Lb (and then the 
linearly related Lx) more than Lfir. Finally, we observe 
how, even taking into account the smallest radius of 6 
Kpc, it is not possible to move the theoretical point near 
the linear relation holding for spirals. 

5.2.2. The masses of stars and gas 

We also investigated in Fig. 1101 what happens if we forget 
about the clues coming from dynamical simulations on the 
masses of stars and gas and we arbitrarily start varying 



12 



G. Galletta et al.: ISM tracers 



the amounts of stars or gas, keeping all fixed. Filled di- 
amonds represent the shift of the model of lowest radius 
if we are changing the mass of stars inside R ga i , going in 
fraction from /» = 0.2 to /» = 1.0, with respect to the to- 
tal amount of stars in the dynamical model. The effect is 
simply to move the point along a line about parallel to the 
linear relation. A smaller amount of stars imply directly a 
lower luminosity Lb (and therefore a lower Lx), but also 
a lower Lpm, because the weaker radiation field makes 
dust cooler and shifts the peak of dust emission to wave- 
lengths longer than 100/ito, with the result of a smaller 
Lfir- Finally, with open circles we show in Fig. 1101 five 
models obtained at fixed amount of stars and at varying 
the mass of gas (and therefore of dust) from fd = 0.2 to 
fd = 1.0, in fraction respect to the total amount of gas in 
the dynamical model. The effect of this huge increase of 
the mass of diffuse gas and dust (in the original model at 
Rgai = 6Kpc only 0.03% of the gas is inside R ga i) is to shift 
the models straight toward the linear relation. It can be 
explained in the following way. Increasing the amount of 
diffuse gas/dust (with all the parameters fixed and the star 
formation exhausted) implies more absorption of the stel- 
lar radiation and therefore a smaller Lb (and Lx)- On the 
other side, Lpm remains almost unchanged or becomes 
smaller. The reason is that the strongly increased mass of 
dust makes the average stellar radiation field weaker and 
therefore the increased emission of dust (due to the big- 
ger mass) peaks at wavelengths longer than 100/im, leav- 
ing Lfir almost unchanged. Even if in this way we can 
shift the model toward the linear relation, the situation is 
physically unrealistic, requiring a huge amount of gas/dust 
concentrated in the centre of an early- type galaxy with ex- 
hausted star formation, which is not commonly observed 
and also not predicted by dynamical models. 



5.2.3. The scale radii 



Further geometrical parameters that must be examined 
are the scale radii r l c of the King's laws - eqn. ljT2l - 
that describe the distribution of the stellar and gaseous 
components. The averaged profiles showed in Fig. and 
used for the models of Figs. and ^| are both charac- 
terized by r* ~ 0.5, allowing for a concentrated amount 
of stars and gas in the inner regions. Keeping all the 
other parameters fixed, we investigated what happens if 
we allow for a uniform distribution of one or both the 
physical components. Three cases have been considered: 
a uniform distribution of gas keeping fixed the stellar one 
(jsas QQ^stars ^ q.5), a uniform distribution of stars 



keeping fixed the gaseous one {r s c 



0.5) 



and, finally, a uniform distribution of both the compo- 

oc). The results are shown in 



nents 



oo,r 



gas 



Fig. ^] for two radii of the galaxy model, R ga i = 6 Kpc 
and Rgai = 20 Kpc, respectively. The three different dis- 
tributions give a similar result: a weaker Lpm, shifting 
the point to the left, and a slightly higher Lx- 




Fig. 11. Models of the more massive galaxy of 1.6 TO 11 M Q 
baryonic mass in the Lpjn vs. Lx diagram for various core 
radii rf ars and r^ as , plotted as filled circles. The thin ar- 
rows indicate the original model, while the thick arrows 
indicate the shift of the model for various core radii. We 
represented: E galaxies (empty diamonds), E/S0 galax- 
ies (pluses), £0 galaxies (crosses), SOs galaxies (empty 
squares), Sa — Sb (empty triangles). 



This can be explained in the following way: for r s J' ars 
and rf as both ~ 0.5, the diffuse ISM and the stars are 
both concentrated in the inner region of the galaxy with 
a density of stars/gas of many order of magnitude big- 
ger than the outer regions. This is the best condition to 
produce high Lfir, because we have that the regions of 
higher density of dust are the same in which there is also 
the higher average radiation field heating dust. The spa- 
tial distribution of the ISM favors the interaction with the 
stellar radiation. When we destroy this coupling between 
stellar emission and density of gas, as we do allowing for 
a uniform distribution of gas or stars or both, the emis- 
sion in the Lfir. becomes weaker. The weakening of the 
dusty emission is stronger for the bigger radius of 20 Kpc 
because in all the three cases one or both the components 
are distributed over a huge galactic volume and we have 
low density of gas eventually coupled with weak radiation 
field. For the 6 Kpc model, even if the coupling in the 
central regions is destroyed, the galaxy is small enough to 
keep a good level of Lpm, even when the matter is equally 
distributed across all the galaxy volume. 

5.2.4. The star formation history 

Last and main point to be examined is how varying the 
star formation history affects the position of the galaxies 
into the Lpm vs Lx plot. In Fies. llOl andll lithe galaxies of 
different morphological type form a sequence that, going 
from systems in which the star formation got exhausted 
long ago to systems in which star formation is still ac- 
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Fig. 12. L_fy,r vs L x for galaxy models of spi rals and star- 
burster s. Five theoretical models, taken from lPiovan et cd\ 
(2006b), are represented: three spirals (M100, M51 and 
NGC6946) and two starbursters (Arp220 and M82). 
Different morphological types are plotted as in Fig. ^2 
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Fig. 13. Different adopted star formation histories for 
the model of L6 ■ W ^Mp, baryonic mass. The SFHs of 
IPiovan et all l|2006bh have been re-calculated for the same 
mass of the dynamical model. The upper panels repre- 
sent two spiral- like SFHs, while bottom panels show two 
starburst-like SFHs, peaked in the last million years. 



tive, moves toward the linear relation and suggests the key 
role played by the star formation. First of all we calculate 
the LprR and Lx obtai ned by the SEDs and the models 
by IPiovan et all l|2006tJ ) of real galaxies of the local uni- 
verse: three spiral galaxies (M100, M51 and 7VGC6946) 
and two starburst galaxies (Arp22Q and M82). The key 
point is that the SFHs of these galaxies allow us to cover 
a good number of different star formation histories. All 
these SFHs, unlike the ones of the ellipticals obtained by 
dynamical simulations, never end and in the case of the 
two starbursters a strong burst of star formation is added 
in the last millions of years. A huge amount of Lfir comes 
therefore from the young and deeply obscured region of 
star formation and not only from the diffuse component. 

The results, presented in Fig. ^] show that the three 
models of spirals stay near the linear relation, while the 
two starbursters stay below the line, with the model of 
Arp220, powered by a huge burst of star formation falling- 
well below the linear relation. The stronger is the emis- 
sion coming from the regions of star formation and the 
bigger is the shift toward higher L fir and lower Lx (due 
to the lower Lb)- The results obtained from the mod- 
els are quite similar to the observational data: for Ml 00 
we get (Lfir, Lx) = (10.28,7.29) with the observations 
giving (10.37,7.01), for Arp220 we have (L FIR ,L X ) = 
(11.92,7.17) compared with (11.99,7.60) and for M82 
we get (L FIR ,L X ) = (10.15,6.45) against (9.79,6.31). 
However, these galaxy models, even if they well represent 
real galaxies, differ in many parameters from the early- 
type galaxy model of 1.6 ■ 10 11 Af Q , like geometry and 
mass. These parameters, together with the SFH, obviously 
concur to determine the position of the models into the 
Lfir vs Lx plot. To isolate the effect of the SFH, we 



first re-calculated the SFHs of the above five theoretical 
models, rescaled to the mass of 1.6 • 10 1:l Mq of the early- 
type galaxy model. In Fig. ^5] we can see four of the five 
SFH obtained. Second, we fixed all the geometrical pa- 
rameters to the same values used for the average model 
of 1.6 • 10 11 M Q early type galaxy. The additional parame- 
ters, that is the escaping time of young stars from parental 
molecular clouds, the library of SEDs of young dusty re- 
gions and the mass of gas in the diffuse a nd molecular 
compon ent, are fixed to the values used in IPiovan et all 
(2006b) for spirals and starbursters as appropriate. 

In Fig. E| we finally show the results obtained as a 
function of the SFH of the galaxy of 1.6 • 10 11 M Q , keep- 
ing all the other parameters fixed. It is interesting to ob- 
serve that since now the star formation never ends and the 
galactic wind is not included, the classical semi-analytical 
chemical evolution can be much more safely coupled to 
the spectro-photometric code. The effect of varying the 
SFH at fixed mass is to enhance the Lfir, keeping al- 
most fixed the Lx and shifting the points toward the 
linear relation at higher infrared luminosities. This is ul- 
timately due to the strong and efficient reprocessing of 
the light coming from very young stars, occurring into 
the dusty star-forming regions. As a consequence of this, 
models with starburst-like SFHs shift, as expected, toward 
higher Lfir luminosity than models with spiral- like SFH, 
because of the stronger star formation and therefore emis- 
sion coming from young dusty regions. This can be also 
understood if we look in detail at the relative contribution 
to L fir. coming from the regions of star formation (let us 
define it Jmcs) and represent it as usual in \og(LFi R / Lq) 
and from the diffuse interstellar medium (Jism)- We get 
the following values: {Jsfr = 10.15, Jism — 10.48), 
(fsFR = 10.02, fisM = 10.45), (f S FR = 9.98 Jism = 
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Fig. 14. Lfir vs Lx for galaxy models of fixed 1.6 • 
lO n M0 baryonic mass for various SFHs (large circles 
and squares) or at fixed SFH for various masses from 

I. 6 • 1O 8 M to 1.6 • 10 n Af© (large filled diamonds). Large 
circles and squares represent respectively models with 
spiral-like SFH and starburst-like SFH. Large filled circles 
and squares are obtained using eqn. while large open 
circles and squares stand upon eqn. ^| Different morpho- 
logical types are plotted as in Fig. ^2 

10.37) for the three models with spiral-like SFHs, while 
we have (f S FR = 10-97, f IS u = 10.47) and (f SFR = 

II. 93, fisM = 10.56) for the models with starburst-like 
SFHs. The stronger is the contribution from star forming 
regions, the higher is Lfir keeping Lb (and Lx) almost 
unchanged. Models slightly dominated by the ISM contri- 
bution, but with a significant contribution coming from 
obscured newly born stars are more suitable to agree with 
the linear relation of spirals. 

It is worth noticing that in Fig. ^] we show both the 
results obtained applying the early-type linear relation be- 
tween L x and L B - eqn. (fTUfl - and the late-type one - eqn. 
£IJ. Since, however, the SFHs used (see Fig. 1131) are typical 
of late type galaxies (or starbursters), it's more physically 
sounded to apply eqn. (J3J to obtain the Lx luminosity. As 
last point we calculated also a sequence of models in which 
one of the SFHs of the spirals has been chosen (namely 
the one of NGC6946) with all the parameters fixed and 
only the mass is varied. As we see from Fig El the effect 
of varying the mass is to shift the object in diagonal al- 
most along the relation. This is simply explained by the 
smaller amounts of stars/gas emitting radiation. 

6. Conclusions 

We have been able to describe the relations existing in a 
galaxy between the various tracers of the ISM and to fix 
the coefficients of the relations existing between FIR, B 
and X-ray luminosity, both for early-type and late-type 
galaxies. 



The large set of data we used allowed us to redefine 
more clearly the relation existing between the CO and the 
10 Mm fluxes. We found that the relation, first obtained 
bv lBreeman et all l|l992f) for early type galaxies, is valid 
also for late type galaxies. In these galaxies, the X-ray flux 
appears linked also to B and FIR emissions. 

The only relation lacking from observations, i.e. the 
one between Lx and Lfir has been studied by the use 
of the most recent chemo-dynamical models coupled with 
dusty evolutionary population synthesis. 

The calculated luminosities of the models seem to con- 
firm our hypothesis about a connection between the ex- 
haustion of the star formation and the " migration" of the 
early type galaxies above the linear relation in the Lx 
vs Lfir plot. In the frame of our assumptions, we may 
therefore conclude that the prediction of our dusty chemo- 
dynamical models of galaxy evolution is consistent with 
the observed lack of a direct relation between Lx and 
Lfir for early type galaxies and is due to the different 
mechanisms of production of FIR light in galaxies where 
the active star formation is no longer active. In most of 
our early-type galaxies the mechanism of IR emission is 
no more strictly related to the ongoing star formation and 
to the reprocessing of the radiation in the dense regions 
where new stars are born. The FIR emission comes there- 
fore most likely from circumstellar dusty shells around 
AGB stars and from an interstellar diffuse medium due 
to the outflow of dusty gas from AGB and RGB stars. 

Finally, we can summarize that: (i) the SFH of the 
galaxies seems therefore to have the stronger effect on the 
position of early- type galaxies in the Lx vs Lfjr plot; 
(ii) other parameters, like the radius of the galaxy and 
the scale radii of stars and gas, play a secondary role, 
even if they can significantly contribute to the scatter of 
the models in the region above the linear relation; (iii) the 
mass is the main parameter explaining the scatter of the 
points along the linear relation. 
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